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Abstract 

An open issue in turbulence theory is related to the determination of the exact evolution equation 
for the probability density associated to the relevant (stochastic) fluid fields. Such an equation 
in the usual approaches to turbulence reproduces, at most in an approximate sense, the correct 
fluid equations. In this paper we present a statistical model which applies to an incompressible, 
resistive and quasi-neutral magnetofluid. The approach is based on the formulation of an inverse 
kinetic theory (IKT) for the full set of MHD equations appropriate for an incompressible, viscous, 
quasi-neutral, isentropic, isothermal and resistive magnetofluid. Basic feature of the new approach 
is that it relies on first principles - including in particular the exact validity of the fluid equations 
- and thus permits the determination of the correct evolution equation for the probability density. 
Specific application of the theory here considered concerns the case of statistically homogeneous 
and stationary MHD turbulence. 

PACS numbers: 47.10.ad,47.27.-i,05.20.Dd 
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I. INTRODUCTION 



The goal of this paper is to investigate a well-known theoretical issue of fluid dynamics. 
This is concerned with the rigorous formulation of the so-called turbulence problem, i.e., the 
description of the dynamics of a classical fluid in the presence of stochastic sources. We refer 
here, in particular, to the statistical description of turbulence for incompressible fluids based 
on its phase-space approach, namely achieved by representing the turbulent fluid means of 
a suitable phase-space probability density function (pdf). In such a context, an open issue 
is related to the determination of the exact evolution equation for the probability density 
associated to the relevant (stochastic) fluid fields in an incompressible turbulent fluid. It 
is well known, in fact, that in customary approaches to turbulence (see for example Monin 
and Yaglom 1| 1975 and Pope, 2000 [2]) this equation reproduces only in an approximate 
sense the correct fluid equations. In this paper we present a statistical theory which applies 
to an incompressible, resistive and quasi-neutral magnetofluid. 

The purpose of this paper is to propose a solution based on the formulation of an inverse 
kinetic theory (IKT) for an incompressible, viscous, quasi-neutral, isentropic and resistive 
magnetofluid. Basic feature of the theory is that it permits the determination of the exact 
evolution equation for the probability density. 



A. Deterministic and stochastic descriptions of fluids 

In fluid dynamics the state of the fluid is assumed to be prescribed by an appropriate set 
of suitably smooth functions {Z} = {Zi, i = l,n} denoted as fluid fields. These are required 
to be real functions which are at least continuous in all points of a closed set Q x I, closure 
of the extended configuration domain fix I. In the remainder we shall require that: 

1. Q (configuration domain) is a bounded subset of the Euclidean space on R^; 

2. I is identified, when appropriate, either with a bounded time interval, i.e., I=]tQ, ti[ C 
R, or with the same time axis (M); 

3. in the open set Qxl the functions {Z} , are assumed to be solutions of an appropriate 
closed set of PDE's, denoted as fluid equations; 
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4. by assumption, these equations together with appropriate initial and boundary condi- 
tions are required to define a well-posed problem with unique strong solution defined 
everywhere in Qxl. 

Depending on the definition of initial conditions (A), boundary conditions (B), possible 
volume forces (C), as well as, more generally, also the functional form assumed for the fiuid 
equations (D), two types of descriptions of the fiuid, deterministic and stochastic, can in 
principle be distinguished, in which the fiuid fields are treated respectively as deterministic 
or stochastic functions (see definitions in Appendix). The two descriptions are obtained, 
respectively, when all conditions (A), (B), (C) and (D) are deterministic, or at least one of 
them [A,B, C or D] is prescribed in a stochastic sense (see below). For definiteness, in the 
following (except for Sec. 7) we shall formally rule out case D, imposing that the functional 
form of the fiuid equations remains the same in both cases. 

B. Two types of turbulence - Homogeneous and stationary turbulence 

The stochastic description of fiuids is usually denoted as turbulence In this respect, 
however, two different types of phenomena should be distinguished, i.e., 

• physical turbulence 

• numerical turbulence. 

The first one arises due to "physical effects", namely the specific definition adopted for 
the fiuid equations and for the related initial and boundary conditions (these definitions 
may be considered, in a sense, as physical characteristics of a fiuid, i.e., as related to the 
properties of the fiuid). Thus, physical turbulence may appear only as a result of the 
stochastic "sources" indicated above. The first two (A and B) are due to possible stochastic 
initial and boundary conditions. In fact, the initial and/or boundary state of the fiuid may 
not be known deterministically, as implied by measurement errors, or by the action of some 
type of initial and/or boundary volume/surface stochastic forcing. Finally (case D), it is 
obvious that the fiuid equations for the stochastic fiuids fields [stochastic fluid equations) 
may, in principle, differ in other ways (besides C), from the corresponding deterministic 
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fluid equations. Their precise prescription is therefore essential for the definition of specific 
turbulence models. 

Instead, numerical turbulence occurs only in numerical simulations and is due specifi- 
cally to approximations involved in the numerical solution method (for the fluid equations). 
Hence, it is necessarily related to the speciflc algorithms adopted for the construction of the 
numerical solutions. 

Both types of turbulence may produce, in principle, "real" effects to be observed in 
numerical experiments. However, it should be stressed that, in practice, the distinction 
between the two phenomena may be difficult, or even impossible, to ascertain in numerical 
simulations. Nevertheless, from theoretical perspective it is obvious that the distinction 
between the two phenomena should be made. 

This paper deals with exact properties and implications of the fluid equations (and the 
related initial and boundary conditions) , disregarding completely issue related to the possible 
choice of the approximate solution methods adopted for their numerical solution. Hence, 
we will face only the flrst type of phenomena. For this reason in the following the term 
"turbulence" will be meant, in a proper sense, only in the physical sense here specifled. 

Turbulence is therefore characterized by the presence of stochastic fluid fields Zi G {Z} , 
for i = l,n (with {Z} deflning the state of the fluid), which are assumed to depend on 
suitable hidden variables a — {ai,i — l,k} e Q K'^, with A; > 1, which are assumed 
independent of {r,t). Hidden variables by assumption cannot be known deterministically (in 
other words, in the context of classical mechanics, they are not "observables"), and hence 
are necessarily stochastic. This means that by assumption: 1) the hidden variables a are 
characterized by a probability density g, denoted as reduced stochastic probability density 
function (rs-pdf) defined on Va , with Va a non-empty subset of C M.'^ and k > 1; 2) that 
the fiuid fields {Z} will generally be considered as stochastic functions (see Appendix). 

Critical issues in turbulence theory are in principle related both to the definition of the 
hidden variables a as well as the determination of the corresponding rs-pdf g. Both depend 
necessarily on the specific assumptions concerning the source of stochasticity. Therefore, 
it is obvious that for both of them the definition is non-unique. In particular, the precise 
definition of g may depend on the specific source (of stochasticity). This is obvious - for 
example - when its origin is either provided by boundary conditions {i.e., surface forces) 
or respectively by volume forces: both sources, in fact, can be prescribed in principle in 
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arbitrary ways. It follows that, in general, g should be considered of the form 



g = gir,t,a), (1) 

with (r,t) G Qxl. This corresponds to the so-called (statistically) nonhomogenous and 
nonstationary turbulence. Particular solutions are, however, provided by the requirements 
that there results identically either g = g{t,a) [homogeneous turbulence), or g = g{r,a) 
[stationary turbulence). 

In subsequent sections (Sec. 2-6) we intend to present a mathematical formulation of 
turbulence in magnetofluids which pertains specifically to the validity of both assumptions, 
i.e., 

9 = 9{c^)- (2) 

This defines the so-called homogenous and stationary turbulence. Possible generalizations of 
the theory to nonhomogeneous and nonstationary turbulence are discussed in Sec. 7. For this 
purpose, all the fluid fields {Z} will be generally assumed stochastic {i.e., non-deterministic), 
leaving unspecified the definition of the hidden variables and - apart for the previous require- 
ment - the form of the rs-pdf. 



C. Turbulence in magnetofluids 

The phenomenon of turbulence in magnetofluids {MHD turbulence) is nowadays playing 
a major role in plasma physics and fluid dynamics research. We refer here in particular to 
incompressible magneto- and conducting fluids for which the fluid description is appropriate 
(for such systems typically a statistical description based on the microscopic dynamics is 
not possible). The widespread picture of MHD-turbulence phenomena occurring in these 
fluids consists of an ensemble of finite-amplitude waves with random phase. Examples of 
;urbulence models of this type are well known. For example, Iroshnikov [sl and Kraichnan 
611 independently assumed that MHD turbulence occurs as a result of nonlinear interactions 
between Alfven wave packets P, |8|. However, there is an increasing evidence that this 
picture is an oversimplification. In particular, MHD turbulence may include fluctuations 
whose phase-coherence characteristics are incompatible with wave-like properties. The latter 
are the so-called coherent structures, like clumps, holes, current filaments, shocks, magnetic 
islands, vortices, convective cells, zonal flows, streamers, etc. 



In fluid turbulence the signature of the presence of coherent structures is provided by the 
existence of non-Gaussian features in the probability density. This is usually identified with 
the velocity-difference probability density function, traditionally adopted for the description 
of homogeneous turbulence. In the past several statistical models have been proposed for 
its determination, which include the mapping-closure method, the test-function method 
and the field-theoretical approach [ll, 12, 13|. Nevertheless, despite the progress achieved 



in modelling key features of the basic phenomenology, still missing is a consistent, theory- 
based, statistical description of MHD turbulence. Clearly, such formulation - if achievable at 
all - should rely exclusively on a rigorous, deductive formulation of the turbulence-modified 
fluid equations following from the fluid equations. 

D. Goals of the investigation 

Based on a recently proposed inverse kinetic theory (IKT) for classical and quantum fluids 
(Ellero and Tessarotto, 2004-2007 QQ, 16| 



equations 



^ V^), and in particular its formulation for MHD 

isl . Q , a phase-space statistical model of turbulence is proposed. The present 



approach, which represents a generalization of the one recently developed for incompressible 
neutral fluids j^, [20I, allows to determine exactly the stochastic evolution equation of the 
local position-velocity joint probability density function {local pdf). The local pdf, rather than 
the velocity-increment pdf traditionally adopted in turbulence theory, is in fact found to be 
meaningful for the dynamics of the fluid. Indeed, we intend to prove that the stochastic fluid 
fields which characterize the turbulent magnetofluid are uniquely determined by means of 
suitable velocity-moments of the local pdf. Nevertheless, as previously pointed out 
IKT's are intrinsically non-unique in character, in particular because they may be also 
formulated in such a way to apply only to a suitable subset of the fluid equations. In this 
respect two types of phase-space approaches can in principle be developed, namely either: 

• a Complete IKT, yielding the full set of stochastic fluid equations which advance in 
time the stochastic fluid fields {Z} ; 

• a Reduced IKT, which applies only to a suitable subset of fluid equations, to be iden- 
tified with the set fluid equations advancing in time only {{Z)} . 
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In particular the Complete IKT can be non-uniquely prescribed (see THM's 1 and 2) so 
that: 



1. the time-evolution of the stochastic fluid fields {Z} describing the state of a turbu- 
lent magnetofluid is uniquely determined by means of a suitable phase-space classical 
dynamical system {the MHD- dynamical system)] 

2. the MHD-dynamical system is uniquely associated to an appropriate inverse kinetic 
equation (IKE) which advances in time an appropriate phase- space probability density 
/ (Pdf); 

3. the inverse kinetic equation (IKE) is such that its velocity moments coincide with 
the complete set of stochastic (MHD) fluid equations. In particular, the IKE's which 
advance in time the stochastic-averaged pdf (/) as well its stochastic fluctuation 5f = 
f — if) , determine respectively the fluid equations for stochastic-averaged fluid fields 
{{Z)} and the stochastic fluctuations {SZ} = {Z — (Z)} ; 

4. under appropriate initial and smoothness conditions, the strict positivity of / is assured 
for arbitrary strictly positive initial pdf's; 

5. under the same assumptions for (/) , exhibits an irreversible time-evolution (see the 
Corollary of THM.2) 

Instead, as previously pointed out [sl, [lol, the Reduced IKT (see THM.3) can be realized 
so that: 

1. the time evolution of the stochastic averaged fluid fields is similarly uniquely deter- 
mined by means of a suitable dynamical system and an appropriate inverse kinetic 
equation advancing in time the stochastic-averaged pdf (/); 

2. the strict positivity of (/) is assured; 

3. the inverse kinetic equation for (/) is Markovian. 



Key feature of present theory is that - unlike several other approaches to be found in 



the 



literature (for a review see for example, Monin and Yaglom, 1975 and Pope, 2000 Q) - 
the relationship between fluid fields and the pdf is exact (and not just only approximate). In 



particular, the evolution of the stochastic-averaged fluid flelds is determined via a suitable 
nonlinear transformation (here generated by an appropriate dynamical system). This is 
achieved by the construction of a Vlasov-type kinetic equation which advances in time the 
local pdf and - as a consequence - also the stochastic fluid flelds of the magnetofluid. 



II. MOTIVATIONS: IKT APPROACH FOR MAGNETOFLUIDS 



A fundamental aspect of fluid dynamics is the construction of phase-space approaches for 
realistic fluids. Indeed, phase-space techniques are well known both in classical and quantum 
fluid dynamics. In fact, generally the fluid equations represent a mixture of hyperbolic and 
elliptic PDE's, which are extremely hard to study both analytically and numerically. This 
has motivated in the past efforts to replace them with other equations, possibly simpler to 
solve or mathematically more elegant. In this connection a particular viewpoint - which 
applies in principle both to classical and quantum fluids - is represented by the class of so- 
called inverse problems, involving the search of a so-called inverse kinetic theory (IKT) able 
to yield the complete set of fluid equations for the fluid flelds, via a suitable correspondence 
principle. As a consequence the fluid equations are identifled with appropriate moment 
equations constructed in terms of the relevant kinetic equation. This raises the interesting 
question whether the theory can be generalized to arbitrary classical magnetofluids. The 
issue is relevant at least for the following reasons: a) the proliferation of numerical algorithms 
in MHD which reproduce the behavior of incompressible fluids only in an asymptotic sense; 
b) the possible veriflcation of conjectures involving the validity of appropriate equations of 
state for the fluid pressure; c) the ongoing quest for efficient numerical solution methods 
for the MHD equations, with particular reference to Lagrangian methods. Finally, another 
important motivation is the possibility of achieving an exact solution to this problem, based 
on inverse kinetic theory (IKT) (see Tessarotto et al, 2008 [l8|). In fact an aspect of fluid 
dynamics is represented by the class of so-called inverse problems, involving the search of 
IKT's able to yield identically a prescribed set of fluid equations. A few examples of IKT's 



which hold both for c 
et al, 2000-2007 14, 



assical and quantum fluids have been recently investigated (Tessarotto 



15 



17 



2l|). 



In particular, among the possible IKT's, special have IKT's which pertain to incompress- 
ible fluids described in terms of strong solutions of the corresponding fluid equations. In 
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such a case IKT's can be defined in such a way that the kinetic equation which advances 
in time the kinetic distribution function generates a suitable classical dynamical system. 
Phase-space approaches of this type have been aheady achieved for several types of fluids 
(Ellero and Tessarotto, 2000-2008 [is], Q, 17, 21]). In particular, the same approach can 



be extended in principle also to incompressible magnetofluids. Without loss of generality, 
in the following we consider the case of the MHD equations describing an incompressible, 
viscous, quasi-neutral, isentropic, isothermal and resistive magnetofiuid. Starting point is 
the assumption that there exists a suitable phase-space classical dynamical system, to be 
denoted as MHD- dynamical system, 

Xo x(t) = Tt^to^o, (3) 
defined in such a way that its evolution operator Tf^to uniquely advances in time the fluid 



fields 



171]. The dynamical system (and the operator Tt^to) is assumed to be generated by a 



suitably smooth vector field X(x, t; Z), solution of the initial- value problem 



d 



X = X(x, t; Z), 



where X(x, t; Z) is required to be - in an appropriate way - functionally dependent on the 
set of fluid fields {Z} which characterize the fluid. Here x = (ri,V]^) G F and respectively 
ri and Vi denote an appropriate " state- vector" and the corresponding "configuration" and 
"velocity" vectors. In particular, it is assumed that x spans a phase-space of dimension 2n 
(r C R^"), where by assumption n > 3. Therefore, introducing the corresponding phase- 
space probability density function (pdf) /(x,t;Z), in F it fulfills necessarily the integral 
Liouville equation 

J(t;Z)/(x(t),t;Z) = /(x„t„,Z), (5) 
to be viewed in the following as a Lagrangian inverse kinetic equation (Lagrangian IKE), 



with J{t; Z) 



dx{t) 



denoting the Jacobian of the map ([3]). In particular, let us assume 



dx.o 

for definiteness that ([3]) is a diffeomorphism at least of class C^^''(r x / x /), with / C M 
denoting an appropriate finite time interval. Then, if the initial pdf f{'Xo,to, Z) is suitably 
smooth it follows that the pdf /(x,t; Z) satisfies also the differential Liouville equation 

L(Z)/(x,t;Z) = 0, (6) 



L{Z) denoting the corresponding streaming operator 

^(^)--|-+|^-{X(x,t;Z).}. (7) 

This equation [and the equivalent Lagrangian equation ([5])] will be considered in the fol- 
lowing as as an Eulerian inverse kinetic equation (Eulerian IKE) for a suitable set of fluid 
equations. In other words, the arbitrariness of the dynamical system ([3]) [i.e., of X(x, t; Z)], 
will be used to seek a representation of the vector field X(x, t; Z) such that that suitable the 
velocity moments of f{'x,t; Z) can be identified with the relevant fluid fields characterizing 
a prescribed fluid. 

Let us consider, for definiteness, the case of an incompressible, viscous, quasi-neutral, 
isentropic, isothermal and resistive magnetofluid subject to the condition of isentropic flow. 
By assumption, the relevant fluids 

{Z} = {p = p,>0,Y,p>0,B,St}, (8) 

i.e., respectively the mass density, the fluid velocity, the fluid pressure, the magnetic field 
and the thermodynamic entropy, are assumed to be defined in the whole existence domain 
Q X I. In particular if Q denotes an open connected subset of R^, its closure Q by definition 
is the set where the mass density is a constant p = po > 0. We shall assume that the 
fluid fields are continuous in i7 x J, satisfy a suitable set of deterministic MHD equations 
in the open set Qxl and, moreover, fulfill appropriate initial and (Dirichlet) boundary 
conditions respectively at t = to and on the boundary of the configuration domain SQ. 
These requirements are represented respectively by the equations 

P = po, 

v-v = o, 

V-B = 

NV = 0, 

< V ■ A^V = 0, (9) 

NbB = 0, 
f 5t = 0, 
Z{T,to)=Zoir), 
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Here the notation is standard. Thus, and A''^ are the nonhnear Navier-Stokes and 
Faraday- Neumann differential operator 

d 1 
N\ = —V + V • VV + — [Vn - fl - z/V^V (10) 

dt Po 

and 

AT^B =— B + V- VB-h--^V2B. (11) 

ot Ana 

Here the vector fields f and h, to be denoted as volume forces, read respectively 

f .p„g+_LB.VB-vg), (12) 
h = -B ■ VV. (13) 

where the mass density po, the kinematic viscosity u and the conductivity cr are all real 
positive constants to be considered non-stochastic. Furthermore, in the first equation, pog 
denotes the gravitational force density and the remaining terms the Lorentz force density. 
In addition, the first three equation in denote respectively the so-called incompressibility, 
isochoricity and divergence-free (for B) conditions. The remaining ones include, instead, the 
Navier-Stokes equation, the energy equation, the Faraday-Neumann equation for B and the 
isentropic entropy condition (with St denoting the thermodynamic entropy). Finally, the 
last two equations denote respectively the initial and Dirichlet boundary conditions. Thus, 
by taking the divergence of the N-S equation, it follows the Poisson equation for the fluid 
pressure p which reads 

V'j9 = -poV • (V ■ VV) + V ■ f, (14) 

with p to be assumed non negative and bounded in Q x I. Finally, in the following and 
consistent with Eqs.Q the thermodynamic entropy St will be assumed as a constant in the 
whole extended configuration domain Q x I. 



III. IKT FOR THE DETERMINISTIC MHD EQUATIONS 



A prerequisite for the subsequent formulation of a turbulence theory based on the phase- 
space approach here adopted [3|] (see subsequent Sections 4-6) is the development of a Com- 
plete IKT for the deterministic fluid equations [see Eqs.(l9])]. As recently pointed out in 
Ref . [l8 1 . this goal can be achieved by introducing suitable generalization of the of the theory 
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earlier developed for the incompressible Navier-Stokes equations [15|, ll6|, ll7|, l2l|]. For defi- 
niteness, let us introduce the notations ri= (r, s) , vi= (v, y) and X = {vi, Fi(x, t; f,Z)}, 
where the vectors r and vi span, respectively, the whole configuration domain of the fiuid 
{Q) and the 3-dimensional velocity space (M'^). Moreover s eQ, y G are two additional 
real vector variables, with s denoting in particular an ignorable configuration-space vector 
[both for the fiuid fields and the kinetic distribution function /(x,t;Z)]. The streaming 
operator L{Z) in this case reads 

where {Z} are the fiuid fields and the vector field 

Fi(x, t; /, Z) = {F(x, t; /, Z), Y(x, t; /, Z)} (16) 

can be interpreted as an effective mean field force acting on a particle state x = (r, s, v, y). 
In the following we intend to prove that, at least in a suitable finite time-interval /, the fiuid 
fields Z{r,t) can be uniquely identified with the velocity moments J (iviG'(x,i(;)/(x,t; Z), 
where /(x,t; Z) is properly defined kinetic distribution function and G(x,t) appropriate 
weight functions. More precisely, there results: 

1 = /dvdy/(x,t; Z), 
V(r,t) =/rfvrfyv/(x,t; Z), 
p,{r,t)=pofdYdylu'f{^,t;Z), (17) 
B(r,t) = /rfvrfyy/(x,t; Z), 

Sr{t) = S{f{t)), 

{correspondence principle). In particular V, p and B are identified respectively with the 
moments G'(x,t) = \',pou'^/3, y, where u = v — V. In addition, if the same distribution 
function f(t) = /(x,t; Z) is strictly positive in the whole set Fx/, S{f{t)) is the statistical 
Boltzmann-Shannon entropy functional 

Sifit)) = -J cix/(x,t; Z) ln/(x,t; Z). (18) 
r 

Consistent with Eqs.Q we shall impose that S{fit)) exists for all t E I, and that there 
results identically (Vt G /) the conservation law 

S{f{t)) = const, (19) 

12 



{constant H-theorem) . To reach the proof, let us first show that, by suitable definition of the 
"force" fields F(x, t) and Y(x, t), a particular solution (which defines a kinetic equilibrium 



35|) of the IKE is delivered by the Maxwellian distribution: 



Here we denote 



= (22) 

'^th,T 

where u = v — V, w = y — B, f^^^ = 2pi/po and f^^j- is an 'a priori'' arbitrary strictly 
positive constant. Furthermore 

pi(r,t) =po(t)+p(r,t) (23) 

is the kinetic pressure. In these definitions, Po(^) (fo be denoted as pseudo-pressure) is an 
arbitrary strictly positive and suitably smooth function defined in /, while the mass m > 
is an arbitrary real constant. The following theorem can immediately be proven: 

Theorem 1 - Complete IKT for the deterministic MHD equations Let us assume 
that in the existence domain Q x I the fluid fields ^ are strong solutions of the MHD fluid 
equations Furthermore, let us identify the vector fields F(x, t;/, Z) and Y(x, t;/, Z) 
with 

F(x,t;/M,^) =Fo + Fi, (24) 



Y(x,t; /m, Z)= u • VB - h--^V'B + wV • B, (25) 

Ana 

where Fq and Fi read respectively 

Fo(x,t; /m, Z) = -i+lu- VV+lvV ■ u+z/V^V, (26) 

Po ^ ^ 

F.(xA/„,Z) = ^^« + fv.np,{|-|}. (27) 

while 
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and ^Pi is uniquely prescribed consistent with the energy equation, namely 
D d 



d 1 
—V^I2 + V • Vl^V2 V ■ f - z/V-V^V 

Ot Po 



(29) 



Finally, let us assume that the initial kinetic probability density f{to) = f{^,to, Z) belongs to 
the functional class {f{to)} , defined so that f{to) fulfils at time t = t^ the moment equations 
[T7\ ), while the initial fluid fields Z{r,to)=Zo{r) are suitably prescribed. 
It follows that: 

1) the local Maxwellian distribution ^2D\j is a particular solution of IKE [Eq.^]\ 

2) this is a solution of the same equation if an only if the fluid fields {Z} satisfy the fluid 
equations 

3) the velocity-moment equations obtained by taking the weighted velocity integrals of 
Eq.^ with the weights 

G{x, t) = Po, V, PomVS, yV/S (30) 

deliver identically the same fluid equations 

4) the initial kinetic probability density f{to) = f{^,to, Z) is assumed to satisfy the prin- 
ciple of entropy maximization (PEM) [34], which requires that 

6Sifito)) = 0; (31) 

5) ^t,to € /, the statistical entropy S{f{t)) satisfies the constraint 

S{f{to)) = S{f{t)) (32) 

{constant-H theorem). 
PROOF 

First, we notice that if the fluid equations (Q are satisfied identically in f2 x /, the proof 
that ( l20l) is a particular solution of the IKE [Eq.(l6])] follows by direct differentiation. The 
converse implication, i.e., the proof that if ( |20l) is a solution of Eq.(l6]) then the fluid equations 
(Q are satisfied identically in Q x I, follows by evaluating in particular the (vi—) velocity- 
space moments of Eq.Q for the weights fl30|) G = po (v — V)^ /3 and w'^/3. These moments 
deliver respectively the two isochoricity conditions 

V • V = 0, (33) 
V ■ B = 0. (34) 
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The proof of Proposition 3) follows by noting that the extremal solution of the varia- 
tional principle fl3T|) coincides necessarily with the Maxwellian distribution /jvf(x, to; Z) [see 
Eg. (1201) ]. Finally, to satisfy Proposition 4) it is sufficient to notice that the pseudo-pressure 
Po{t) [which enters the definition of the kinetic pressure (1251) ] can always we defined in such 
a way to satisfy, at least in a finite time interval /, the constraint (I32p . Q.E.D. 

In analogy to the case of isothermal fiuids the present theorem can be generalized in a 
straightforward way to a suitably smooth non-Maxwellian initial distribution function [l^. 



IV. STOCHASTIC MHD FLUID EQUATIONS 

Let us now construct the stochastic MHD fiuid equations appropriate in the case of homo- 
geneous and stationary turbulence. In principle this requires the adoption of an appropriate 
turbulence (or stochastic) model. This concerns the introduction of suitable hypotheses, i.e., 

1. the definition of an appropriate set of hidden variables a; 

2. the definition of the associated rs-pdf g; 

3. the assumption that at least some of the fluid fields, together with their initial and 
boundary values on 6Q, and/or the volume forces f and h depend in a suitable way 
from the hidden variables a; 

4. the prescription of the form of the stochastic fluid equations. 

Its should be stressed that in principle the definition of the a's - and consequently of 
g and of the dependencies on a to be suitably prescribed in the fiuid equations - remains 
completely arbitrary. Regarding item 4), in particular, it is obvious that 'a priori ' also 
the stochastic fiuid equations might differ from the deterministic equations iQ. Thus, for 
example, the condition of incompressibly of the fluid p = Po might be violated in a turbulent 
fluid due to the possible presence of stochastic mass density fluctuations (see Sec. 7). 

For deflniteness, let us assume that the fluid flelds are stochastic, i.e., they depend 
smoothly on the hidden variables a, while the stochastic fluid equations retain the same 
functional form of the corresponding deterministic equations ([9]). This implies that the fluid 
remains, in particular, incompressible, viscous, quasi-neutral, isentropic, isothermal and re- 
sistive magnetofluid, which - in particular - implies necessarily 6p=0 and 6St = 0. Hence 
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the meaningful fluid fields are only represented by the subset 



{Z{a)} = {V(r,t, a),p{r,t, a) > 0, B(r,t, a)} . 



(35) 



Together with their initial and boundary values on 6Q (defined respectively by the vector 
fields Zo{r,a) and Z{r,t,a)\g^ 

we shall require that they all admit the stochastic decompositions: 



Zw{r,t,a)\g^), and the volume forces f(r,t, a), h(r,t,a), 



Z(r± a) 



Zo(r; a) 



Z^{r,t,a) 



f(r,t,a) 



{Z{r,t, a)) + 5Z{r,t, a), 
{Zo{r,a)) +5Zo{r,a), 
{Z^{r,t,a)) + SZ^{r,t,a) 
(f(r,t,a)) +5f(r,t,a). 



(36) 
(37) 
(38) 
(39) 



Here the stochastic-averaging operator () is defined according to Eg. (1771) . while (Z) , {Zo) , 
{Zw) , (f) and 6Z, 6Zo, 6Zw and 5f are denoted respectively as stochastic averages and 
stochastic fluctuations. 

In the remainder the precise definition of the operator () [i.e., the specification of g{a) 
and of the hidden variables a] is not needed. The only requirement to be introduced is 
that - consistent with the assumption ([2]) - it commutes with all the differential operators 
appearing in the MHD equations, in particular, ^, V and V^. In such a case the fluids 
equations for the averaged fields {{Z{a))}, are immediately found to be 



(p) = Po, 
V • (V) = 0, 

V-(B) = 0, 

(N) (V) + {6N6V) = 0, 
(V) • [{N) (V) + {5N5\)] = 0, (40) 
{Nb) (B) + {5Nb6B) = 0, 

I (St) = 0, 
{Z{r,to;a)) = {Zo{r;a)) , 
^ {Zir,t;a))\g^= {Z^{r,t; a))]^^ , 
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which are denoted as the stochastic-averaged MHD equations for {{Z)}. Analogous equa- 
tions hold for the stochastic fluctuations {5Z{a)} : 

V ■ 5V = 0, 

V ■ 5B = 0, 

(A^) 5V + 5N5Y~ {SNS\) = 0, 
< (V) ■ [{N) 6Y + 6N6V- {6N6Y)] = 0, (41) 
(Nb) 6B + 6Nb6B- {6Nb6B) = 0, 
6Z{r,to; a)= 6Zo{r; a), 
6Z{r,t;a)\s^ = 6Z^{r,t;a)\s^ 

{stochastic-fluctuating MHD equations ). Here the notation is standard [isl . Finally, 
(N) , 6N and respectively (Nb) , SNb are the nonlinear operators 

(^) (V) = ^ (^) + (^) ■ V (V) + 1 [V {p) - (f )] - uV' (V) , (42) 



and 



SN6V = 6y ■ V5V + — [V6p - 6{] - uV^SV, (43) 

Po 

{Nb) (B) = I (B) + (V) • V (B) - (h) -^V^ (B) , (44) 

(JiVs^B = 5V • V5B - 5h - ^V^^B. (45) 

Ana 



V. COMPLETE IKT FORMULATION FOR THE STOCHASTIC MHD EQUA- 
TIONS 

As indicated above (see also related discussion in Ref. [sl) in principle the problem of 
the formulation of a stochastic IKT can be posed either for the full set of fluid equations or 
only for the stochastic-averaged equations fHOl) (see also Sec. 6). 

In the first case let us assume, for greater generality, that the pdf can depend on the hidden 
variables both implicitly, via the fiuid fields, and possibly also explicitly, i.e., it is stochastic 
local pdf of the form /(x,t; Z{a), a). The formulation of the IKT in this case follows directly 
from THM.l by imposing that /(x,t; Z{a),a) satisfies the stochastic Eulerian IKE 

L{Z{a))f{^,t; Z{a), a) = 0. (46) 
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Here, the streaming operator L{Z{a)) is again defined by Eq. (fT5!l . while by assumption the 
vector fields F and Y depend on a only implicitly via the fiuid fields {Z{a)}. Next, let us 
introduce for the pdf and the streaming operator the stochastic decompositions 

fit; a) =/(x,t; Z{a), a)= (/(x,t; a)) + 5/(x,t; Z{a), a), (47) 

L{Z{a)) = {L{Z{a))) + 6L{Z{a)), (48) 

where () is again the stochastic-averaging operator (1771) . Requiring that () commutes also 
with the differential operators d/d\ and d/dy, Eq. fl76|) yields the two stochastic kinetic equa- 
tions advancing in time respectively (/) = (/(x,t; Z{a),a)) and 6f = 5/(x,t; Z{a),a), i.e. 

{L{Z)) if) + {6L{Z)6f) = 0, (49) 

(L(Z)) 6f + 5L{Z)5f - {5L{Z)5f) = 0. (50) 

Thanks to THM.l, Eq.( l46l) [or the equivalent Eqs. fl49l) and (1501) ] provides, in principle, the 
exact (and unique) evolution of the stochastic pdf /(x,t; Z(a), a). Furthermore, thanks 
to the correspondence principle [given by Eqs. ffTTl) ]. also of stochastic fluid flelds {Z} are 
uniquely determined j^. 

In this Section we intend to formulate a Complete IKT, yielding the full set of stochastic 
fluid equations advancing in time the stochastic fluid flelds {Z} , here represented by the 
stochastic MHD problem deflned above. The basic result has the flavor of: 

Theorem 2 - Complete IKT for the stochastic MHD equations Requiring that 
the fluid fields {Z{a)} = {V,p > 0, B} are stochastic [in the sense the of the definition given 
above], and the correspondence principle 

1 = J dYdyf{yi,t;Z{a),a), 
V(r,t; a) = f dvdy\f{yi,t;Z{a),a), 
< pi(r,t) = po/dvdy|u2/(x,t;Z(a),a), (51) 
B(r,t;a) = f dYdyy{'x,t; Z{a),a), 
STit;a) = Sifit;a)), 

holds identically in Tx I (where T is the closure of T = Qx M^^) then provided the mean-field 
force 

Fi(x, t; f, Z{a), a) = {F(x, t; f, Z{a), «), Y(x, t; f, Z{a), a)} (52) 
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is still defined by Eqs. (fl^- ( f^) (with fluid fields to be considered stochastic functions) it 
follows that: 

1) the local distribution fM{'^,t; Z{a)), i.e., the Maxwellian distribution ^2D) with stochas- 
tic fluid fields [or, equivalent, the contributions {fM{^,t]Z{a))) and 6fM{:K,t;Z{a))] is a 
particular solution of the inverse kinetic equation [respectively of Eqs. (|7^-( r5^) / if only 
if the fluid fields {{Z{a))} and {6Z{a)} satisfy respectively the stochastic fluid equations 
(g^ and (gip; 

2) fM{^,t; Z{a)) maximizes the Boltzmann-Shannon entropy S{f{t; a)) in the functional 
class {/(x, t; {Z{a)) ,a)} defined so that the pdf satisfies solely the constraints (5l\) . where 
the stochastic- averaged fluid fields {Z{a)} are considered prescribed; 

3) the velocity-moment equations obtained by taking the weighted velocity integrals of 
Eq.^[4^ [or, equivalent, Eqs. (f7^-(f7^/ '"'^^^ the weights l[3U\) deliver identically the same 
stochastic fluid equations and ( fT?] ). 

PROOF 

The proof is an immediate consequence of THM.l and follows by invoking the stochastic 
decompositions ( H7I) and ( 148|) . together with the assumption ([21) for the rs-pdf g, which 
assures the equivalence of the two stochastic kinetic equations P9|) - fj49|) with Eq.Q. In 
particular, it is immediate by direct evaluation to verify that stochastic fluid equations (HOj) 
and fHT]) are recovered from Eqs. (1491) - fHOl) by taking their moments fl30|) . Q.E.D. 

Main consequences of THM.2 are that: 

1. the inverse kinetic equation, either represented in the Eulerian or Lagrangian forms. In 
particular, the Eulerian equation is provided by Eq. (H^ [or by the equivalent Eqs. (H^ 
and (!50l) ]. while the Lagrangian equation, in view of Eq.([5]), is again of the form ([5]), 
namely 

J (t; Z{a), a)/(x(t, a),t; Z{a), a) = /(x^, to, Z{a), a) (53) 

where Xq x(t, a) is the stochastic MHD-dynamical system determined by the 
(stochastic) phase-space Lagrangian equation (jl]) with vector field X(x, t; Z(q;)), 



J{t; Z{a), a) 



Ox(t) 
9xo 



denoting the corresponding Jacobian. 



2. the two phase-space representations, provided respectively by the Eulerian and La- 
grangian IKE's, are equivalent. In fact, by construction, both equations satisfy iden- 
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tically the stochastic MHD equations, including the appropriate initial and boundary 
conditions imposed on the fluid flelds; 

3. the time-evolution of the stochastic pdf / is, in turn, uniquely determined by the 
stochastic dynamical system deflned above [see, in fact, the Lagrangian inverse kinetic 
equation Eq. (!53!) ]. 

4. The stochastic-averaged pdf (/) can be proven to be strictly positive. This result 
is important to assure both that (/) is truly a probability density and that its time 
evolution is (possibly) irreversible. This result can be shown to be a consequence of 
THM.l. 



VI. MARKOVIAN REDUCED IKT FOR THE STOCHASTIC MHD EQUATIONS 



In the traditional approach to turbulence theory main emphasis is related to the construc- 
tion of the statistical equation advancing in time the stochastic-averaged pdf, rather than 
on the complete pdf. An example is provided by stochastic models - based on Markovian 
Fokker-Planck models of small-scale fluid turbulence recently investigated in the literature 
by several authors (including: Naert et al., 1997 22|; Friedrich and Peinke et ai, 1999 |23| : 
Luck et ai, 1999.24; Cleve et al, 2000 [251; Ragwitz and Kantz, 2001 [261]; Renner et ai, 



2001, 2002 
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28|; Hosokawa, 2002 |29|). 



An interesting issue is whether the possible Markovian character of the statistical equation 
advancing in time the stochastic-averaged pdf, in the present case to be identifled with the 
local pdf (/) , can be established based on an inverse kinetic theory approach. 

As previously pointed out Q, the stochastic IKE provided by Eq. fj46|) . or the equiv- 
alent equations ( l49l) and ( !50l) . are formally similar to the Vlasov equation ari sing in the 



kinetic theory of quasi-linear and strong turbulence for Vlasov- Poisson plasmas |30|, l3l|, l32| 
and related renormalized kinetic theory 33|], which are known to lead generally to a non- 



Markovian kinetic equation for (/) alone. Nevertheless, the stochastic- averaged kinetic 
equation [i.e., Eq. fj49|) ] is known to be amenable, under suitable assumptions, to an ap- 
proximate Fokker-Planck kinetic equation advancing in time (/) alone. This is achieved by 
formally constructing a perturbative solution of the equation ( l50l) for the stochastic pertur- 
bation 6f. To obtain a convergent perturbative theory, however, this usually requires the 
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adoption of a suitable renormalization scheme in order to obtain a consistent kinetic equa- 
tion for (/) . This raises the interesting question whether in the framework of IKT there 
exist possible alternatives based on the construction of exact Markovian IKT formulations. 

A solution to this problem is provided by the construction of an IKT only for the 
stochastic-averaged equations fHUl) . i.e., to be achieved by means of a Reduced IKT j^. 
Here we intend to show that this leads to an inverse kinetic equation of the form ([6]), 

L((Z(a)))/(x,t;(Z(a))) = 0, (54) 



which is manifestly Markovian [3|, ll9| and satisfies also a constant H-theorem (i.e., it con- 
serves the Boltzmann-Shannon entropy). Here / = /(x, t; (Z)) is the Eulerian local pdf for 
the stochastic INSE problem, which advances in time the stochastic-averaged fluid fields 
{{Z{a))} and L{{Z{a))) is the corresponding streaming operator, to be defined in terms of 
the vector field 

Fi(x, t; /, {Z{a))) = {F(x, t; f, (Z(«))), Y(x, t; /, (Z(a)))} , (55) 
with Fi to be interpreted again as a mean field force acting on a particle with state x. 



Besides the requirement of validity of the fluid equations fHUl) |l6l ]. let us impose that vector 
X has the form X = {v, y,F,Y}, where F(x,t; /, (Z(a))) and Y(x,t; /, (Z'(a))) are both 
assumed generally functionally dependent on the local pdf. Again, by appropriate choice of 
the mean field forces F, Y, the moment equations can be proven to satisfy identically the 
fluid equations [Eqs. fl40l) ]. as well the appropriate initial and boundary conditions. Requiring 
that / is a strictly positive let us require that there results identically in F x / : 

l = JdYdyf{^,t- (Z)), 
(V(r,t;a)) = /rfvrfyv/(x,t; (Z)), 
pi(r,t) = p J dvdyl {nf /(x, t; {Z)) (56) 
(B(r,t;a)) = /rfvrfyy/(x,t; (Z)), 
{STit;a)) = Sifi^,t;{Z))). 
Here ^(/(x, t; (Z) ,a)) is the (Shannon) entropy integral 

5(/(x, t; {Z) ,«)) = - ^ c?x/(x, t; {Z{a))) In /(x, t; (Z(a))), (57) 

u = V — V(r, t] a) is the relative velocity, while 

Pi{T,t)=Po{t) + {p{r,t-a)) (58) 
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is the kinetic pressure and Po(t) a smooth real function (pseudo-pressure) to be suitably 
defined (see below). Again the form of the local pdf / = /(x, t; {ZY) can be chosen in such 
a way to satisfy the principle of entropy maximization (PEM) j34], i.e., imposing, at the 
initial time t = to, the variational equation SS{f) = 0, with 6'^S{f) < 0, while requiring that 
/ belongs to the functional class {/(x, t; {Z{a)))} , manifestly prescribed by the constraints 
(!56|) . Hence, PEM implies yields necessarily (for t = to) that the initial pdf must be of the 
form 

/M(x,t; {Z)) = exp {-X^ - Y^} . (59) 

Here, in difference to Eqs. fl21l) and fl22|) . X"^ and now denote 

-^^ - 

= (61) 

^th,T 

where u = v — V, w = y — B, f^^^ = 2pi/po and v^f^ j. is an 'a priorV arbitrary strictly pos- 
itive constant. Eq. fl54|) implies the construction of a suitable classical dynamical system, 
defined by a phase-space map Xq — > ii.{t) = Tt^to^o, where Tj^^ is the evolution opera- 
tor [l6| generated by the stochastic phase-space Lagrangian equation (jll) with vector field 
X(x, t, {Z{a))) and Jacobian J(x(t),t, {Z{a))). 

Therefore, the Eulerian IKE can also be represented in the Lagrangian form [equiv- 
alent to Eq.fl54 



J(x(t),t, (Z(a)))/(x(t),t; {Z{a))) = /(x„,t,; (Z,(a))) = /„(x„; (Z„(a))), (62) 

{Lagrangian IKE), Here, /o(xo; {Zo{a))) is a suitably smooth initial pdf and 

ax(t) 



J{^{t),t;{Z{a))) 



(63) 



<9xo 

is the Jacobian of the map Xq x(t). Then it is immediate to prove that the stochastic 
MHD problem admits a Reduced IKT. For the sake of simplicity here we present the result 
which holds in the case of a Maxwellian (local) pdf of the type (l59ll . 

Theorem 3 - Markovian Reduced IKT for the stochastic-averaged MHD equa- 
tions 
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Let us assume that: Ai) the stochastic MHD problem [J^ j admits a smooth strong solution 
inT X I, with I a finite time interval; A2) the mean-field force 

Fi(x, t; f, {Z{a))) = {F(x, t; f, Y(x, t; f, {Z{a)))} (64) 

is defined by 

F(x,t;/M,(Z(a)))=Fo + Fi, (65) 



Y(x,t; /m, {Zia)))= (u ■ VB) - (h) - (B) + (wV • B) , (66) 

where Fq and Fi read respectively 

Fo(x,t; /m, {Z{a))) = 1 (f ) + 1 (u ■ VV) +1 (VV • u) (V) , (67) 

Po ^ ^ 

Fi(x,t; /m, {Z{a))) = + ^Vlnpi - ^1 , (68) 

^ = I (Poit) + (P)) + (V) -V (Po(t) + (P)) (69) 

is prescribed by the requirement of validity of the stochastic- averaged energy equation \l^J; 
A3) /(x,t; {Z{a))) satisfies suitable initial and boundary condition consistent with the initial- 
boundary value problem ( f^^ (see Ref. A4) the initial local pdf, f{x,to] {Zo{a))), is 
suitably smooth an strictly positive. It follows that: 

1) the velocity-moment equations of IKE [54\ ) evaluated for the weight functions (^(x, t) = 
V, coincide with the stochastic INSE equations) ; 

^) ^/ f{to, is strictly positive so it is f{to, {Z{a))) in the whole domain Fx/; 

3) the Maxwellian local pdf ( f5P|) is a particular solution of IKE if an only if the fluid 
fields {{Z)} satisfy the stochastic-averaged MHD problem (^^; 

4) the pseudo-pressure Po{t) can be uniquely determined in the time interval I in such 
a way that, denoting fit, {Z{a))) = /(x, t; (Z)), in the same time interval I there results 
identically 



pifit, {Z{am = (70) 



(constant H-theorem). 
PROOF 
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The proof is immediate for proposition 1), while 2) follows invoking the Lagrangian IKE 
fl62|) and by noting that by construction J(x(t),t; {Z{a))) > in the whole domain Fx/. 
To prove 3) let us assume that a strong solution of the stochastic INSE problem exists. 
In such a case it is immediate to prove that /A/(x(t),t; (Z)) is a particular solution of the 
inverse kinetic equation (15^ . In fact, substituting (1591) in the inverse kinetic equation (15^ 
it follows: 

L/M(x,t; {Z{a))) = {| (V) +v ■ V (V)| •^/m+ 
|(B).v.V(B)}.ig);„ 

+ (|lnpi+v ■ VlnpA - l]fM- (71) 
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which manifestly implies Eqs.(l40l). Instead, if we assume that in F x /, / = fui'^i t; {Z{a))) 
is a particular solution of the inverse kinetic equation (IKE) ( 15^ . which fulfills identically 
the constraint equation (170]) . thanks to proposition 1) it follows that the fluid fields (V) , (p) 
are solutions of the INSE equations. Finally to prove proposition 4) and 5), let us invoke 
IKE to evaluate the entropy production rate, which reads 



^5(/(t, {Z{am = - I dx/(x,t; {Z)) 



r 



|;-F(x,t;/) + |^-{Y(x,t;/)} 



(72) 



Hence, in a bounded time interval I it always possible to satisfy the constraint placed by 
the constant H-theorem, requiring 

|5(/(t, {Z{am = (73) 



d'^r- 



|(Po(t) + (p)) + (V)-V(Po(t) + (p)) 



0. (74) 



2 Jn Po{t) + {p) 

This delivers an ordinary differential equation for the pseudo-pressure Po{t). Assuming that 
the fluid fields (V) , (p) are suitably smooth, this equation can always be fulfilled at least in 
the case in which the domain of existence / is a finite time interval. Q.E.D. 

Let us briefly analyze the consequences of the theorem, with particular reference to the 
condition of strict positivity of the pdf [proposition 2)] and the constant H-theorem [Eq. (173!) . 
see proposition 5)]. It is well-known that H-theorem assures, for arbitrary (and suitably 
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smooth) initial conditions of the pdf, also the strict positivity of the kinetic distribution 
function. Nevertheless, in view of the validity of proposition 2), the requirement of its 
validity is, in principle, unnecessary. In fact, the H-theorem is only necessary because of the 
requirements posed by the condition of isentropic flow [see Eqs. fl40l) ] and the correspondence 
principle [defined by Eqs.f l56p ]. However, there is (another) important physical consequence 
of the theorem. In fact, if the requirement is posed that the initial pdf /(to, {Z{a))) satisfies 
PEM, i.e., it maximizes S{f{to, {Z{a)))), the constant H-theorem assures that in the the 
domain of existence J, f{t, {Z{a))) must necessarily maximize S{f{t, {Z{a)))). Therefore, if 
for all t & I, fit, {Z{a))) is assumed to belong to the same functional class specified above 
{fit, {Z{a)))} , then necessarily fit, (Z(a))) is a kinetic equilibrium, to be identified with 
the local Maxwellian distribution flMJl) . Similar conclusions apply manifestly to THM.2 too, 
provided the functional class is suitably defined (see THM2.) and the kinetic equilibrium is 
identified with (([20])). 

Another interesting issue is provided by the comparison between the two stochastic IKT- 
approaches discussed above, i.e., the Complete and the Reduced IKT's here developed. 
Key differences between them are manifestly provided both by the form of the inverse ki- 
netic equations and the choice of the initial conditions for the corresponding pdf. This 
is due to the different definition of the mean field force in the two caseSs, given respec- 
tively by Eqs. fl52l) . flMl) and moreover by the different requirements posed in the two case 
by PEM (the principle of entropy maximization; Janyes, 1957 34|). In fact, the initial pdf 
/(to, {Zi<y))), which in both cases is an extremal of the Boltzmann-Shannon entropy, namely 
is such that 6Sifito, {Zia)))) = 0, belongs to different functional classes {/(to, (Zia)))}. 
The precise definition of {/(to, {Zia)))} in the two cases is actually important since it de- 
termines uniquely the initial condition [see related discussion in Refs jSj and |20|]. Provided 
either the total initial fiuid fields Z ir, to] a)=Zoir; a) or only the stochastic averaged ones 
(Z(r,to; «)) = {Zoir; a)) are prescribed 

In both cases [see THM's 1, 2 and 3] the initial kinetic probability density is provided by 
a Maxwellian-type kinetic equilibrium distribution. This is found as the "most-likely" initial 
kinetic probability density at t = to, in agreement with PEM. In particular, in the first case 
the distribution /a/(x, t;Z) [see Eq. (l20|) ] follows by requiring that {/(to, a)} is prescribed 
imposing the constraints defined by Eqs. fll7l) - to be evaluated at t = to - while in the second 
/a/(x, t; {Z)) [see Eq.(l59l)] follows from the constraints (l56l) . This means that in the first 
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case the initial fluid fields the compete set of fluid fields {Z}^^^^ must be considered as 
prescribed, while in the second only the stochastic-averaged fluid fields { (Z) is initially 
known (so that in this case the initial value of the stochastic fluctuations {SZ}^^ remains 
in this case arbitrary). 



VII. GENERALIZATIONS TO NONHOMOGENEOUS AND NONSTATIONARY 
TURBULENCE 

Turbulence models depend both on the definitions of the hidden variables a and of the 
related rs-pdf g (see also previous section). This involves specifically the case in which g 
is taken of the form g{r,t]a), which corresponds to the assumption of nonhomogeneous 
and nonstationary turbulence. It is obvious that, unless specific additional assumptions are 
introduced, such models are non-unique, due to the arbitrariness of their possible choices 
[both for a and g]. Their unique determination, however, may in principle be achieved either 
based on suitable 'ad /ioc' mathematical models or the phenomenology of (magneto-)fluids, 
based either on first principles [i.e., implied by the underlying microscopic molecular dynam- 
ics and statistical mechanics), or the observation of real fluids or obtained from numerical 
experiments. Nevertheless, despite this arbitrariness, in the particular case in which the 
functional form of the fluid equations [i.e., Eqs.dH])] remains unchanged, the generalization 
of the present theory is still possible. This requires solely suitable smoothness assumptions 
to be satisfied by the rs-pdf. On the contrary, a similar conclusion is not obvious if the 
form of the fluid equations [Eqs.([9])] is modified. In fact, in such a case, i.e., for arbitrary 
turbulence- modified fluid equations, the existence of an IKT-approach 'a priori' cannot 
generally be assured. 



VIII. CONCLUDING REMARKS 



A statistical model of MHD turbulence has been pointed out utilizing the inverse kinetic 
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2l|. Basic 



theory (IKT) recently developed for classical and quantum fluids Q, Q, 0, 
feature of the theory is that the IKT has been constructed in such a way to satisfy exactly the 
stochastic MHD equations describing an incompressible, viscous, quasi-neutral, isentropic. 



isothermal and resistive magnetofluid. As consequence an inverse kinetic equation has been 
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determined which uniquely advances in time the local pdf. 

The main result of the paper is represented by THM.2 and 3, yielding IKT- approaches 
to the stochastic MHD equations and particularly by the inverse kinetic equations [i.e., 
Eqs.f l49p . flSU]) and dSlD], which advance in time respectively the stochastic or the stochastic- 
averaged local pdf. The first two equations [i.e., Eqs.f HU]) and fHI]) ]. determine uniquely the 
evolution of the stochastic fluid fields, i.e., both the stochastic averages {{Z{a))} and the 
corresponding stochastic fluctuations {5Z{a)} . Instead, Eq. fl54l) determines only {{Z{a))} 
{i.e., provided {5Z{a)} is known). In both cases the inverse kinetic equations have the prop- 
erty that their velocity-moments (to be deflned in terms of appropriate weighted velocity- 
space integrals) satisfy identically a closure condition. In fact, there exists, by construction, 
a subset of the moment equations which is closed. Such a set is deflned so that it coincides 
with the appropriate stochastic fluid equations [i.e., either Eqs. fHOl) or fHT]) ]. 

The theory displays several interesting new features, with respect to customary phase- 
space statistical approaches to turbulence (see for example . In particular: 

1. it is based on the introduction of a local position- velocity joint probability density 
function (local pdf), rather than the usual velocity-difference pdf (which typically 
affords only an approximate description); 

2. Eulerian and Lagrangian phase-space formulations are manifestly equivalent; 

3. a remarkable, and in a sense, surprising feature of the present theory is that the 
equation advancing in time the stochastic-averaged pdf is can also be set in the form of 
a Markovian kinetic equation This feature is actually "built in" the IKT approach, 
thanks to the deflnition here adopted for the time evolution operator Tt^to which 
advances in time the local pdf (and hence the same fluid flelds). 

In our view the IKT-theory here presented is a useful setting for the investigation of 
theoretical and mathematical aspects of turbulence phenomena which may potentially occur 
in a variety of fluids (both incompressible and compressible) and in particular magnetofluids. 
Several interesting applications and generalizations of the theory are possible, which include 
- for example - the investigation of turbulence in fluid mixtures and in thermo-magnetofluids. 
The proper treatment of such phenomena requires a consistent description of phase-space 
dynamics, to be based on the IKT approach developed in this paper. We remark, in this 
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connection, that in the present approach both the specification of the hidden variables a and 
of the reduced stochastic probabihty density (g) remains in principle completely arbitrary. 
This feature is potentially important since it permits, in principle, the systematic treatment 
of all the possible different sources of stochasticity here pointed out, a problem which remains 
still fundamentally unsolved to date. Finally, the extension of the theory to nonhomogeneous 
and nonstationary turbulence has been pointed out (Sec. 7). Related developments will be 
discussed in greater detail elsewhere. 
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IX. APPENDIX: STOCHASTIC AND DETERMINISTIC FUNCTIONS 

A possible convenient definition of stochastic function A (to be identified with the generic 
fluid field Zi) is the following one: 
Definition - Stochastic function 

A function A is denoted as stochastic on the set fl x I x Va if: 

1) ^4 = A{r, t, a) is assumed to be a real function defined on the set Q x 7 x T4 which 
depends on the real hidden- variables (to be denoted as stochastic parameters) a e 14 C R'^; 

2) the a's are characterized by a probability density g, denoted as reduced stochastic 
probability density function (rs-pdf) defined on , with a non-empty subset of C 
and k > 1. Hence, g is necessarily non-negative {g > 0) [or even strictly positive {g > 0)] 
and such that 




(75) 



da denoting the canonical measure on V^, 
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3) introducing the stochastic averaging operator 

(■) = [ dag; (76) 



where the integration is performed at constant (r, t) G x J, the stochastic function A{r, t, a) 
is assumed to admit everywhere in f2 x / the stochastic average 

(A) = j dagA{T,t,a); (77) 

In contrast, a deterministic function can be defined as: 
Definition - Deterministic function 

A function A{r,t,a) is denoted as deterministic if there results identically on the set 

X Fq, X / 

A={A). (78) 

This requires that that either A is independent of the parameters a or the probability density 
g is a /c-dimensional Dirac delta, i.e., it is of the form g = 5^''\a — «„) = Yl ^i'^i ~ Oioi), 

i=l,k 

with ao a prescribed element of V^. 
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